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*0 . Abstract 
i> . 

We have performed a new model independent amplitude analysis of reaction 
ir p — > 7r 7r + n measured at CERN at 17.2 GeV/c on polarized target. The 
analysis extends the range of dipion mass from 600-880 MeV used in the pre- 

Qh' 

D ' vious study to 580-1080 MeV. Our purpose is to study the role of /o(980) 

interference and to better determine the resonance parameters of <r(750) reso- 

•i— I , 

nance from fits to both S"-wave amplitudes [S'pS and |5| 2 S with recoil nucleon 
transversities "up" and "down", respectively. This work is the first attempt 
to determine resonance parameters of a resonance from Breit-Wigner fits to 
production amplitudes of opposite transversity. We used several fitting pro- 
cedures to understand the resonant structure of the S'-wave amplitudes in 
this mass range. First we performed separate fits to \S\ 2 Y, and |5"| 2 E using a 
Breit-Wigner parametrization for <r(750) and /o(980), and a constant coher- 
ent background. We obtained an excellent fit to the data with a low x 2 P er 
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data point (x 2 /dpt). However the width of a is narrow 110-76 MeV in |S| 2 £ 
while it is broad 217-233 MeV in |5| 2 S. To decide whether the apparent 
dependence of T a on nucleon spin is a real physical effect or an artifact of 
separate fits, we performed simultaneous fits to |5| 2 S and |5| 2 S. First we 
assumed a common single a pole. We found a good fit with m a = 775 ± 17 
MeV, r CT = 147 ±33 MeV but the x 2 /dpt is larger than x 2 /dpt from separate 
fits. We concluded that the results of separate fits may indicate existence 
of two a poles, one with a narrow and the other with a broad width. Self- 
consistency requires that both poles contribute to both S'-wave amplitudes. 
The simultaneous fit to |S*| 2 S and \S\ 2 Ti with two common a poles yields an 
excellent fit with a lower x 2 /dpt than the separate fits. All amplitudes are 
dominated by a a state with m a = 786 ± 24 MeV and T a = 130 ± 47 MeV. 
All amplitudes receive a weaker contribution from a narrower a' state with 
77V = 670 ± 30 MeV and Y c , = 59 ± 58 MeV. We propose to identify a' (670) 
and a (786) with colour-electric and colour-magnetic modes of lowest mass 
scalar gluonium ++ (gg). 
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I. INTRODUCTION 



In 1978, Lutz and Rybicki showed that almost complete amplitude analysis of reac- 
tions ttN — > 7t + tt~N and KN — > K + tt~N is possible from measurements in a single exper- 
iment on a transversely polarized target. More recently it was shown [fJ,[J that amplitude 
analyses of reactions ir~p — > 7i°Tr°n and n~p — > r]7T~p are also possible from measurements on 
transversely polarized target. The work of Lutz and Rybicki opened a whole new approach 
to hadron spectroscopy by enabling to study the production of resonances on the level of 
spin amplitudes rather than spin-averaged cross-sections. 

In amplitude analyses of these reactions resonances are observed in mass distribution 
of recoil nucleon transversity amplitudes \A\ 2 T, and |A| 2 S corresponding to recoil nucleon 
spin "up" or "down" relative to the scattering plane. In traditional hadron spectroscopy, 
the nucleon spin is irrelevant and resonance should appear in both spin amplitudes with the 
same mass, with the same width and with similar mass distributions. The measurements of 
reactions ix~p — > 7i^7i + n, 7i + n — > 7i + 7i~p and K + n — > K + ix~ P at CERN on polarized targets 
seem to challenge these expectations. These CERN experiments reveal shapes of resonant 
mass distributions in amplitudes of opposite transversity that are different and dependent 
significantly on the nucleon spin 



An important advantage of amplitude analyses is that new resonances can be detected 
on the level of spin amplitude that are not visible in the spin-averaged measurements on 
unpolarized targets. Such is the case of the new scalar resonance cr(750) found in amplitude 
analyses of reactions n~p — > 7i^7i + n at 17.2 GeV/c at low momentum transfers —t = 0.005- 
0.2 (GeV/c) 2 and in ir + n — > 7t + tt~p at 5.98 and 11.85 GeV/c at larger momentum transfers 
—t = 0.2-0.4 (GeV/c) 2 |16H T^|. Using Monte Carlo method of amplitude analysis, it was 
found in Jl7j that in n^p — > 7r~7r + n reaction both solutions for spin "up" S'-wave amplitude 
\S\ 2 H resonate below 880 MeV while both solutions for spin "down" S'-wave amplitude \S\ 2 T, 
appeared nonresonating. This result is in perfect agreement with x 2 minimization method 
of amplitude analysis used in [|],[5| (see Fig. 1 and 2 of fll8"| ). Detailed quantitative Breit- 
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Wigner fits to |5| 2 S in |Tj| showed the importance of coherent background for determination 



of the resonance parameters of a(750). The best fit in JTS] gives m a = 753 ± 19 MeV and 
T a = 108 ± 53 MeV. The excellent agreement of these values for m a and T a with the 



Ellis-Lanik theorem [19[] strongly suggests that cr(750) is the lowest mass scalar gluonium 

Tornqvist expressed concern f20| , |2~T|| that the low mass and the narrow width of cx(750) 
is due to the neglect of interference of <r(750) with fo (980) in our Breit- Wigner fits which 
initially extended only to 880 MeV in dipion mass using the data set of Ref. [5]. The 
apparent nonresonating behaviour of the amplitude \S\ 2 H below 880 MeV was also puzzling. 
To investigate these questions we have performed a new amplitude analysis ofn~p— > n^n^n 
reaction at 17.2 GeV/c and small momentum transfers —t =0.005-0.20 (GeV/c) 2 extending 
the dipion mass range up to 1080 MeV using the data set of Ref. [7]. 

In this report we present the results of our detailed study of resonant structure of S- 
wave amplitudes \S\ 2 H and \S\ 2 H in the mass region of 580-1080 MeV. This work is the first 
attempt to determine resonance parameters of a resonance from fits to production amplitudes 
of opposite nucleon transversity. All previous determinations of resonance parameters used 
spin-averaged mass distributions or, in the case of Ref. [18], only one spin amplitude (|S| 2 E). 

The new results for \S\ 2 Ti shows a narrow resonance below 880 MeV, followed by an 
enhancement above 900 MeV and a dramatic dip at ~1000 MeV corresponding to /o(980). 
The amplitude \S\ 2 H shows a broader structure with a dip at ~1000 MeV. 

First we performed separate fits to |5| 2 E and \S\ 2 using a coherent sum of a(750) and 
/o(980) Breit- Wigner amplitudes and constant background. We find that a contributes to 
both amplitudes but with different widths. Solutions 1 and 2 of amplitude |S| 2 E have a 
narrow width of 100 and 76 MeV, respectively. Solutions 1 and 2 of amplitude \S\ 2 H have a 
broader width of 217 and 233 MeV. The interference with /o(980) is crucial in reproducing 
the enhancements above 900 MeV and the dip at 1010 MeV. 

The narrow width in the spin "up" amplitude |S| 2 E and the broad width in the spin 
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"down" amplitude |S| 2 £ is inconsistent with the traditional view of a resonace as a particle 
with single mass and width. To see if we can describe the mass spectra of |S| 2 £ and |S| 2 £ 
with a single a with the same resonance parameters, we have performed a simultaneous fit 
to \S\ 2 H and |S| 2 £ hoping that the differences in the resonant shapes can be accounted for 
by differences in the background. We obtained an acceptable fit but the x 2 P er data point 
{x 2 /dpt) was significantly higher than x 2 /dpt in separate fits. We thus obtained a poorer 
fit. 

Our next step was to admit that a narrow a contributes to |S| 2 £ and a broad a con- 
tributes to |S| 2 £. But a self-consistency requires that the narrow o contributes also to |S| 2 £, 
and the broad a contributes also to |S| 2 £. We expected these additional contributions to be 
small in a simultaneous two-pole fit to |S| 2 £ and |S| 2 £. The results of the fit were surpris- 
ing. Our best fit had x 2 /dpt better than the separate fits to |S| 2 £ and |S| 2 £. And the fit 
required in all amplitudes the presence of two a resonances: higher mass and broader a with 
average mass m a = 786 ±24 MeV and IV = 130 ±47, and a lower mass and narrower a 1 with 
average mass of m a i = 670 ± 30 MeV and IV = 59 ± 58 MeV. In all S-wave amplitudes the 
higher mass a dominates the a' contribution. Since this solution gives a very low x 2 /dpt, 
it must be taken seriously. We propose to identify the <r'(670) and <r(786) as color electric 
and color magnetic modes of lowest mass scalar gluonium ++ {gg), respectively, following 
the gluonium classification scheme of Bjorken |22-24[ elaborated in |25| . 

The paper is organized as follows. In Section II we present and discuss the results of our 
new amplitude analysis. In Section III we introduce the parametrization of unnormalized 
S'-wave amplitudes with a single a pole and present the results of separate fits to amplitudes 
\S\ 2 H and |S| 2 £. Simultaneous fits to these S'-wave amplitudes using a single a pole are 
presented in Section IV. In Section V we introduce the parametrization with two a poles 
and discuss the results of simultaneous fits to |S| 2 £ and |S| 2 S using this parametrization. 
In Section VI we discuss the possible constituent nature of the two a states a' (670) and 
cr (786) resulting from the best fit. The paper closes with a summary in Section VII. 
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II. AMPLITUDE ANALYSIS 



The high statistics CERN-Munich measurement of n p — > n 7i + n at 17.2 GeV/c on 
polarized target was reported in four data sets with kinematics given in the following Table: 

Set —t m Am Rcf. 

(GeV/c) 2 (MeV) (MeV) 

1 0.005-0.20 600-900 20 [4,5] 

2 0.01-0.20 580-1780 40 [6] 

3 0.005-0.20 580-1600 20 [7,10] 

4 0.20-1.00 620-1500 40 [8,9] 

In the Table Am is the size of the mass bins. There is only one t-bin covering the whole 
interval of indicated —t. Only report the results for amplitudes for the Set 1. In 
and ]10|] ratios |S|/|S| and S-wave intensity Is are given, so it is possible to reconstruct 



approximately amplitudes |S| 2 E and |S| 2 E for the Sets 2 and 3. These are shown in Fig. 5 
of Ref. [18] and Fig. 2 below, respectively. All amplitude analysis [[|-TD] use y 2 minimization 
method to find solution for amplitudes and their errors. 

For invariant masses below 1000 MeV, the dipion system in reactions ttN — > tt + tt~N is 
produced predominantly in spin states J = (S'-wave) and J = 1 (P-wave). The exper- 
iments on polarized targets then yield 15 spin-density-matrix (SDM) elements, or equiva- 
lently 15 moments, describing the dipion angular distribution |fE|,|l4| . The measured normal- 
ized observables are expressed in terms of two S'-wave and six P-wave normalized nucleon 
transversity amplitudes. In our normalization 

|S| 2 + |S| 2 + |L| 2 + |L| 2 + \U\ 2 + |I7| 2 + |Af + |iV| 2 = 1 (2.1) 

where A = S, L, U.N and A = S, L,U, N are the normalized nucleon transversity amplitudes 
with recoil nucleon transversity "down" and "up" relative to the scattering plane. The S'- 
wave amplitudes are S and S. The P-wave amplitudes L, L have dimeson helicity A = 
while the pairs U, U and N, N are combinations of amplitudes with helicities A = ±1 
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and have opposite t-channel-exchange naturality [PJl4||. The unnatural exchange amplitudes 
receive contributions from 'V and u Ai v exchanges. The natural exchange amplitudes N, 
N are both dominated by " 'A2 exchange. 

Amplitude analysis expresses analytically the eight normalized moduli 

\S\,\S\, \L\,\L\,\U\,\U\,\N\,\N\ (2.2) 

and six independent cosines of relative phases 

cos(7sz / ),cos(7s C /),cos(7 L ( 7 ) (2.3) 



cos(7 5L ), cos{j su ), cos(7 L[/ ) 



in terms of the measured observables In (2.3) cos 7 and cos 7 are cosines of relative 



phases between pairs of amplitudes with opposite transversity (e.g. between S and L, and 
between S and L). 

There are two similar solutions in each (m, t) bin JT1JIH- However in many (m, t) bins the 
solutions are unphysical: either a cosine has magnitude larger than 1 or the two solutions 
for moduli are complex conjugate with a small imaginary part. Unphysical solutions also 
complicate error analysis. Two methods are used to find physical solutions for amplitudes 
and their errors. They are y 2 minimization method and Monte Carlo method. The two 
methods are described in P ]T7j . [l"8| . 



In Ref. [17] we used Monte Carlo method for finding the physical solutions for amplitudes 
and their errors in — > 7r~7r + n using the data Set 1 for dipion masses in the range 600- 
900 MeV. In Ref. [18] we show that the results obtained using Monte Carlo method agree 
with the results using \ 2 minimization method but are smoother and give consistently lower 
X 2 /dpt in Breit-Wigner fits. 

In this work we report results of Monte Carlo amplitude analysis of 7t~p — > 7T~7T + n 
reaction using the data Set 3 in the mass range 580-1080 MeV in order to study the <r(750) — 
/o(980) interference. Above 1000 MeV our analysis is only approximate as we neglect the 



P-wave. However, below 1080 MeV the P-wave contribution is still small and its neglect 
does not affect the structure of S-wave amplitudes. This is confirmed by direct comparison 
with S'-wave amplitudes obtained in |10[ using the \ 2 method with D-wave included above 
980 MeV (see Fig. 2d below). The Monte Carlo method is based on 40,000 selections of 
spin-density matrix elements within their errors in each (m, t) bin. No physical solution was 
found in 3 bins - 650 MeV, 850 MeV and 930 MeV. 

The amplitude analysis is carried out for normalized amplitudes \A\ 2 and \A\ 2 , A = 
S, L, U, N. However, from the spectroscopic point of view the relevant information is con- 
tained in the unnormalized amplitudes |A| 2 E and |A| 2 S where £ = cPa/dmdt is the reaction 
cross-section. It is the unnormalized amplitudes |v4| 2 £ and |^4| 2 £ which represent the direct 
spin-dependent contributions to the mass distribution cPa/dmdt. The unnormalized moduli 
|A| 2 E and |A| 2 S were calculated using the £ = d 2 o jdmdt from Fig. 12 of Ref. [26]. 

The results of our new Monte Carlo amplitude analysis of 7r~p — > -K~-K + n in the mass 
range 580-1080 MeV are shown in Fig. 1. The first thing we observe is that the p° peak 
in da jdmdt is not uniformly reproduced in all P-wave amplitudes. In fact, the amplitudes 
|iV| 2 £ and |P| 2 £ show considerable suppression of p° production. The shapes of mass 
distributions with opposite nucleon transversities show considerable differences in both S- 
and P-wave amplitudes. The amplitudes (S^E, |£| 2 £, |P| 2 S with recoil nucleon transversity 
"down" are smaller and broader than the amplitudes |S| 2 £, |£| 2 £ and |?7| 2 S with recoil 
nucleon transversity "up". The opposite is true for the natural exchange amplitudes |iV| 2 £ 
and |iV| 2 £. 

The S'-wave amplitude |S| 2 S shows a resonant behaviour below 900 MeV in both solu- 
tions. The relative phase 7 Si between S and L amplitudes is near zero in Solution 1 and 
a small constant in Solution 2. Since |£| 2 £ clearly resonates, the amplitude |S| 2 S must 
resonate as well. The amplitudes |S| 2 S and \L\ 2 H are both broader and their relative phase 
7sl is again near zero in Solution 2. This strongly suggests that the amplitude |S| 2 £ also 
resonates. 

Both solutions in both S-wave amplitudes dip at 1010 MeV. This dip is accompanied 
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with a dramatic change of phase 7 5i and 7 SC/ above 1000 MeV. This behaviour is interpreted 
as evidence for a narrow resonance / (980). Its contribution must be taken into account in 
fits to amplitudes \S\ 2 T, and (S'pS. 

In Figure 2 we summarize the results for the S'-wave amplitudes. Fig. 2a and 2b are based 
on Monte Carlo method of amplitude analysis of data Set 1 and 3, respectively. The ampli- 
tude \S\ 2 H clearly resonates below 880 MeV in both solutions but shows an enhancement 
between 880 MeV and 980 MeV which is particularly large on Solution 2. The enhance- 
ment is followed by a rapid drop at 990 MeV with a dip at 1010 MeV. The apparent broad 
structure of |S| 2 E amplitude is due to poor separation of the resonant structure below 880 
MeV and the enhancement between 880 and 980 MeV which are of comparable size. The 
rapid drop of the enhancement at 990 MeV and dip at 1010 MeV are clearly seen even in 
\S\ 2 H amplitude. As we shall see later the enhancement is due to interference of / (980) 
with coherent background. 

Fig. 2c and 2d are based on x 2 minimization methods of amplitude analysis of data Sets 
1 and 3, respectively. Fig. 2c is based on Fig. 10 of Ref. [5] and Fig. 2d is based on Fig. 2 
of Ref. [10]. The results are consistent with the Monte Carlo amplitude analysis. 

The spin-averaged S-wave intensity is defined as Is = (|S| 2 + |<S'| 2 )E. Since there are 
two independent solutions for the amplitudes |<S'(i)| 2 and |S(j)| 2 , i,j = 1,2, we obtain 4 
solutions for I s which we label as follows 

Is(i,j) = (\S(x)\ 2 + \SU)\ 2 )^ (2-4) 
The results are shown in Fig. 3. We see that only solutions is(l, 1) and is(2, 1) are clearly 
resonating while the solutions 15(1,2) and I ${2, 2) do not have a clear resonant structure 
due to enhancement in the range 880-980 MeV. It is also for this reason that it is better to 
study directly the resonant structure of the amplitudes |S| 2 E and |S| 2 £ rather than their 
spin averaged S'-wave intensity. 

III. SEPARATE FITS TO AMPLITUDES \S\ 2 Z AND |S| 2 £ 
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A. Parametrization of amplitudes 



To understand the resonant structure of the S'-wave amplitudes \S\ 2 T> and \S\ 2 H we first 
performed separate fits to these amplitudes using a parametrization which included single 
<t(750) Breit-Wigner pole, /o(980) Breit-Wigner pole and a coherent background. Following 
Ref. [18], we write the unnormalized S'-wave amplitudes in the form 

\S\ 2 i: = q\F\ 2 (3.1) 
\S\ 2 E = q\F\ 2 

where q is the mass dependent phase space factor determined in [18] to be simply the c.m.s. 
momentum of 7r~ in the 7r~7r + rest frame. The amplitudes F and F are unnormalized ampli- 
tudes which are parametrized in terms of resonance contributions and coherent background. 
For amplitude F we write 

F = R a {s, t, m)BW a {m) + R f (s, t, m)BW f {m) + Q{s, t, m) (3.2) 

where BWr is the Breit-Wigner amplitude 

BW R = 2 m f . (3.3) 

where ttir is the resonant mass, R — a,f. In the following / will refer always to /o(980) 
resonance. The mass dependent width T(m) depends on spin J and has a general form 

r = Ts( ±y+iEdS*l (3.4) 

q R Dj(qr) 

In (3.4) q R = q(m = m R ) and Dj are the centrifugal barrier functions of Blatt and Weishopf 



D (qr) = 1.0 (3.5) 

D x {qr) = 1.0 + {qrf 
10 



where r is the interaction radius. We recall that 

Re BW R = ~f )\BW R \ 2 = w R \BW R \ 2 
m R l 

Im BW R = \BW R \ 2 (3.6) 

In (3.2) the term Q(s,t,m) is the coherent nonresonant background. The amplitude F has 
a form similar to (3.2) with obvious replacements R a — > R a , Rf — > Rf and Q — > 

The energy variable s is fixed and will be omitted in the following. Since the experimental 
mass distributions are averaged over a broad t-bin —At = 0.005-0.20 (GeV/c) 2 , we will work 
with amplitudes averaged over the same t interval. The t-averaged amplitudes have the same 
form as (3.2), so we will simply assume that R a , Rf and Q are all t-independent. We will 
also assume that R a , Rf and Q depend only weakly on the dipion mass m and can be taken 
as constants in our fits. It is then convenient to factor out the phase of R a and define 

R a = IR^ = y/Nge* (3.7) 

C = R f /R a = d+ iC 2 

B = Q/R a = B x + iB 2 
Then the unnormalized amplitude (3.2) has the form 

F = yfN s {BW a + CBWf + B}e^ 
and the parametrization of \S\ 2 H reads 

\S\ 2 Z = qN s {[l + 2w a B 1 + 2B 2 ]\BW a \ 2 + 

+Bf + B\ + \C\ + C%]\BW f \ 2 + 

+2[(w a \BW a \ 2 + BJiwfd - C 2 ) + 

+ (\BW a \ 2 + S 2 )(Ci + w f C 2 )]\BWf\ 2 } 
The parametrization of \S\ 2 T, has the same form as (3.9). 
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(3.8) 



(3.9) 



B. Results of separate fits to |S| 2 E and |SfE 



The measured mass distributions \S\ 2 Y1 and \S\ 2 H were augmented at 650 and 850 MeV 



mass bins by data points from Monte Carlo amplitude analysis using the data Set 1 [17 



Independent fits were performed to both solutions of both S-wave amplitudes with mass 
and width of a resonance being free parameters. The fitting was done using the CERN 
optimization program FUMILI ||28j| . Initially the /o(980) was fixed with mass mf — 980 
MeV and width Tf = 48 MeV. An improved fit was obtained when rrtf and V f were allowed 
to be free parameters to be fitted. The best fit was obtained with fitted /o(980) using as 
initial values of parameters B i} C i: i = 1, 2 the output values from the fit with fixed / (980). 
The results are shown in Fig. 4a. The fitted values of m CT , T a ,mf and V f and the associated 
X 2 /dpt (x 2 per data point) for each solution are given in Table I. 

We see from Fig. 4a that the parametrization (3.9) reproduces well the structure of S'- 
wave amplitudes |S| 2 £ and |S| 2 £ in both solutions. The a resonance peaks are clearly visible 
in both solutions to \S\ 2 Y, while the presence of a resonance in |S'| 2 S is required to explain 
the broad structure observed in this amplitude. The interference of /o(980) with coherent 
background and the a resonance is responsible for the enhancements in the mass region 
880-980 MeV, the rapid decrease at 990 MeV and the dip at 1010 MeV in all solutions. 

We see from the Table I that the masses m a are similar in all solutions. We might expect 
the same situation for the width r CT but this is not the case. The width of a from fits to 
\S\ 2 Y> is a narrow 110 or 76 MeV. The width of a from fits to |>S| 2 S is much broader at 217 or 
233 MeV. We obtain a puzzling result that the width of a depends on nucleon transversity. 

It is an accepted view that hadron resonance is a particle that, upon its production in a 
hadronic reaction, freely propagates and decays independently of other particles involved in 
the reaction. Thus its width cannot depend on the spin of recoil nucleon. However, so far this 
assumption has not been tested in experiments as all meson widths were determined from fits 
to spin-averaged mass distributions. This work is the first attempt to determine resonance 
parameters from different spin dependent mass distributions. It is therefore necessary to 
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determine whether the apparent dependence of IV on nucleon transversity is an artifact 
created by separate fits to \S\ 2 H and |S , | 2 5H, or whether it is a genuine property of resonances. 
To this end we turn to simultaneous fits of \S\ 2 Y1 and \S\ 2 T< examining these two possibilities. 
Finally we note that we have performed separate fits to \S\ 2 E and \S\ 2 Y1 obtained by x 2 
minimization method (Fig. 2d). The fits are very similar to results for Monte Carlo method 
(Fig. 4a) but the x 2 jdpt is larger for all solutions. 

IV. SIMULTANEOUS FIT TO \S\ 2 T, AND |5| 2 S WITH A SINGLE a POLE 

Hadron resonance in a production process is represented as a pole in helicity or transver- 
sity amplitudes with mass and width independent of helicities or transversities of external 
particles. In this Section we examine the possibility that the same a resonance contributes 
to both transversity amplitudes F and F (see (3.8)), and that the difference in the shape of 
|S| 2 E and \S\ 2 T> is due to difference in coherent background. Thus we will try to show that 
the apparent dependence of T a on nucleon spin is an artifact of separate fits to |5| 2 E and 
\S\ 2 H distributions. 

We use the parametrization (3.9) for \S\ 2 H and \S\ 2 T, with common m a and T a and use 
the programme FUMILI [^8| to perform a simultaneous fit to data on \S\ 2 H and |>S| 2 X. The 
two solutions for \S\ 2 H are labeled 1 and 2, while the two solutions for \S\ 2 H are labeled 1 
and 2. Thus there are 4 simultaneous fits (1, 1), (1, 2), (2, 1) and (2, 2). Initially the /o(980) 
was fixed with a mass mj — 980 MeV and width = 48 MeV. An improved fit was obtained 
when irtf and V f were allowed to be free parameters. The best fit was obtained with fitted 
/o(980) using as initial values of parameters Bi,Ci, B^C^i = 1,2 the output values from 
the fit with fixed /o(980). The results are shown in Fig. 4b. The fitted values of m a , F a , mf 
and Tf and the associated \ 2 /dpt for each combination of solution are given in Table II. 

We see in Fig. 4b that the simultaneous fits are similar to separate fits in Fig. 4a. Also, 
the cross-fits (12) and (21) are similar to fits (11) and (22), respectively. In Table II we 
notice that all fits have similar values of m a , T a ,rrif and Tf. The average values of mass and 
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width of a resonance are 



m a = 775 ± 17 MeV , T CT = 147 ± 33 MeV (4.1) 



This compares with m a = 753 ± 19 MeV and T a = 108 ± 53 MeV obtained in [To) from 
fits only to |5| 2 S without the interference with /o(980). The simultaneous fit of \S\ 2 J2 and 
\S\ 2 Y> and the interference with / (980) thus require a slightly higher mass and a broader 
width of a state. 

We now compare the x 2 /dpt for the separate and simultaneous fits: 



Fit 


(11) 


(12) 


(21) 


(22) 


Average 


Separate 


0.250 


0.177 


0.267 


0.194 


0.222 


Simultaneous 


0.386 


0.292 


0.370 


0.272 


0.330 



In this table the x 2 /dpt is the average value of x 2 /dpt for solutions i and j, i,j = 1,2 
calculated by the program FUMILI. We see that \ 2 1 dpt for simultaneous fits is systematically 
larger than for separate fits. The average x 2 /dpt over all combinations is 0.222 for separate 
fits and 0.330 for simultaneous fits, or 50% larger. To reject the hypothesis that Y a depends 
on nucleon spin we would need a simultaneous fit with average \ 2 '/dpt over all combinations 
smaller than 0.222. Since this is not the case, we will explore in the next section the 
possibility that r CT depends on nucleon spin in another simultaneous fit to |5| 2 S and \S\ 2 H 
At this point we comment on our use of x 2 ■ Let A" be the number of observations, M 
the number of free parameters and v = N — M the number of degree of freedom. Then \ 2 
per degree of freedom, x 2 /doj = xV^ i s used to decide the goodness of the fit |29]. The 
fit is said to be good if x 2 /dof<l. In comparing various fits it is more appropriate to use 
X 2 per data point, \ 2 /dpt = x 2 '/N, calculated by the program FUMILI. x 2 /dpt reflects the 
absolute value of x 2 obtained and fits with lower x 2 '/dpt are said to be better than fits with 
higher x" 2 ''/dpt. Notice that x 2 /dof = (x 2 /dpt)(N/u). We have verified that all our fits are 
good with x 2 /dof < 1. 
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V. SIMULTANEOUS FIT TO |Sf £ AND |S| 2 £ WITH TWO a POLES 



A. Two-pole parametrization 

Let us accept that the amplitude \S\ 2 H has a a pole with a narrower width and that the 
amplitude \S\ 2 H has a a pole with a broader width. We can express the helicity amplitude 



in terms transversity amplitudes. Using the notation [IS 



\S\ 2 Y, = q\F\ 2 = q\\F - iF x \ 2 (5.1) 



2' 



1 



\S\ z Y, = q\F\ 2 = q-\F + iF l \ 2 
the helicity nonflip and flip amplitudes can be written 

Fo = 75 (F+T) (5 ' 2) 

The helicity amplitudes acquire two a poles from the transversity amplitudes F and F. Let 
us tentatively label m w and T w the mass and width of the a pole in amplitude |5"| 2 S, and m a 
and T a the mass and width of the a pole in the amplitude [S^E. Inverting the relationship 
(5.2) we get 

F = -^(F + tF 1 ) (5.3) 
F = ±={F -iF l ) 

The self consistency now requires that the transversity amplitudes F and F also possess the 
two sigma poles a and a. We could imagine that while a dominates jS^S, the a presents 
a small contribution. Similarly a could dominate |5*| 2 S with a being a small contribution. 
If this were the case then the apparent dependence of a width on nucleon spin could be 
explained in terms of two sigma poles. Such explanation would still be in accord with the 
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standard picture of hadron resonances provided that the resulting \ 2 /dpt would be smaller 
than the x 2 '/dpt for separate fits. 

Consequently we will write the amplitudes F and F in the following form: 

F = R W BW W + R a BW a + RfBWf + Q (5.4) 

F = R a BW a + R W BW W + RfBW f + Q 

where BW^, BW a , BWf are Breit-Wigner amplitudes for a, a and / resonances and Q, Q 
are background terms. Since we expect the terms BW& and BW a to be dominant in \S\ 2 H 
and l-S'pE, respectively, we can factor out the residue terms R& and R a 

F = \R W \{BW W + DBW a + CBWf + B}e** (5.5) 

F = \R a \{BW a + DBWf + CBWf + B}e i4> 

In (5.5) D, D, C, C, B, B are complex numbers which we shall write in the general form 

D = D 1 + W 2 , ...,B = B 1 + iB 2 (5.6) 

We will again assume t-averaged amplitudes (5.4). This makes the parameters (5.6) t- 
independent and we will assume also their mass independence. The parametrization for 
\S\ 2 H can be written 

[S^Y, = Y id + Y new (5.7) 

where the Y u does not contain the terms involving the parameters D\ and D 2 and has the 
same form as (3.9). With a notation 

H a = \BW«\ 2 , H a = \BW a \ 2 , H f = \BW f \ 2 (5.8) 

we get 

You = qN s [(l + Iw^ + 2B 2 )H a + b\ + b\ 
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{C\ + c\)H f + 2{(w a H r7 + B 1 )(w f C 1 - C 2 ) 



+ (ff (T + J B 2 )(C 1 + W/ C 2 )}# / ] 



(5.9) 



Y 



new 



qN s [2{w a H a + B 1 ) + (C lW/ - C 2 )i/ / }{ J D 1 «; CT - D 2 }H a + 



(5.10) 



+2{(# (T + B 2 ) + (d + ty / C 2 )^ / }{Di + D 2 w a }H a 



+{D 2 i + ^2}^] 



In (5.9) and (5.10) the coefficients w a , w a and Wf are defined in (3.6) for R = a, a and /. 
The normalization term N$ = \Rw\ 2 - 

The parametrization for \S\ 2 H is similar to (5.9) and (5.10) with some changes: the free 
parameters are without bars, and the role H a , W a is interchanged with H a and w a . 



We have fitted the parametrizations (5.9) and (5.10) for \S\ 2 T, and \S\ 2 T, simultaneously 
to the data on \S\ 2 H and \S\ 2 H. Since there are two solutions for each distributions \S\ 2 H 
and |/S r | 2 5H, we obtain 4 independent fits which we label (11), (12), (21) and (22). 

We started with the expectation that the poles a in \S\ 2 H and a in \S\ 2 T, will not change 
much, and that a and a will be only small contributions to \S\ 2 H and l-S'pE, respectively. 
The actual results were different and surprising. 

We run the optimization program with well over 20 initial conditions. In each case the 
fits to combinations (21) and (22) converged to the same solutions for resonance parameters 
m w , IV, m a and IV. Since these solutions are very similar, we call them Solution A. The fits 
to combination (11) resulted in a different solution for the resonance parameters which we 
call Solution B. The fits to combination (12) produced two solutions very close to solutions 
A and B. In each combination (ij) the solution for resonance parameters produced always 
the same x 2 /dpt. 



B. Results 
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Even with additional runs we failed to find Solution A in the combination (11) in free 
unconstrained fits. However when we run a constrained fit with m w , IV, m a and r CT from 
the Solution A from fits to (12), (21) and (22) combinations, we obtained excellent fits 
with x 2 /dpt nearly equal to x 2 /dpt for Solution B from combination (11). The resonance 
parameters from the fit to (12) produced the closest x 2 /dpt. We have therefore accepted 
this constrained fit as the Solution A for the (11) combination. 

Similarly we were unable to find Solution B in unconstrained fits to combinations (21) 
and (22). When we run a constrained fit with m w , IV, m a and T a from the Solution B from 
(11) combination, we obtained acceptable fits with somewhat larger x 2 /dpt in comparison 
with x 2 /dpt from Solution A. On the basis of x 2 /dpt alone the Solution A is much preferable. 

The fits from the Solution A and Solution B are shown in Fig. 4c and Fig. 4d, respectively. 
The figures show two fits to each solution of amplitudes \S\ 2 H and (S'pS. One fit comes from 
the diagonal combinations (11) or (22) (solid lines) while the other is from the off-diagonal 
combinations (12) or (21) (dashed lines). In Fig. 4c and Fig. 4d we see that the two a poles 
do not result in two peaks in the \S\ 2 T, spectrum. However, there is a hump at 670 MeV in 
\S\ 2 Z in Solution A (Fig. 4c) and a small peak at 710 MeV in |S| 2 £ in Solution B (Fig. 4d). 

Our initial expectation was that the a and a resonances will dominate the \S\ 2 H and 
\S\ 2 T, amplitudes, respectively. Two features of a and a poles have emerged in all fits. First, 
a has always a higher mass and a broader width than a. Second and surprising feature is 
that a is always the dominant contribution in all fits to both amplitudes \S\ 2 H and \S\ 2 H 
while a presents a weaker contribution to both amplitudes. Thus our initial expectation 
that the amplitudes and \S\ 2 T, are dominated by a narrow a and a broad a states is 

not validated. Instead, both mass spectra of \S\ 2 Y1 and \S\ 2 T, can be described in terms of 
dominant a and subdominant a. This parametrization with two new a poles leads to x 2 /dpt 
that is lower than x 2 /dpt for separate fits (see below). Thus the hypothesis that the width 
of a depends on nucleon spin must be abandoned as an artifact of separate fits. 

Since W and a poles are no longer associated with the nucleon transversities of the 
amplitudes \S\ 2 T, and l-S'pE, we will relabel these states in a more conventional way. We 
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will make the following change of notation 



a — > cr, ttv -> m CT , I> -> T a 

a -> cr', m CT -> rrv, T CT — >• IV (5.11) 

The results for ttv T a (old r?v, IV) and r?v, IV (old m CT , IV) are shown in the Table III. 
The average values of resonance parameters over the combinations give for Solution A 

m a = 786 ± 24 MeV, T a = 130 ± 47 MeV (5.12) 

77V = 670 ± 30 MeV, IV = 59 ± 58 MeV 

Solution B is similar but the differences between the masses of a and a' are smaller. The 
average values of resonance parameters over the combinations are 

m a = 768 ± 24 MeV, T a = 128 ± 53 MeV (5.13) 

77V = 711 ± 44 MeV, IV = 70 ± 78 MeV 

The Table III also shows the x^/dpt for each combination (ij) which is the average 
of x 2 /dpt for \S\ 2 H and \S\ 2 H solutions. The lowest x 2 /dpt is always obtained for the 
(22) combination. The largest x 2 1 Apt is always obtained for the (11) combination. It is 
worthwhile to summarize the x 2 /dpt of various fits 



Fit 


(11) 


(12) 


(21) 


(22) 


Average 


Separate 


0.250 


0.177 


0.267 


0.194 


0.222 


Simult. (lcr) 


0.386 


0.292 


0.370 


0.272 


0.330 


Simult. (2a) 












- Solution A 


0.274 


0.168 


0.172 


0.110 


0.181 


- Solution B 


0.268 


0.168 


0.246 


0.180 


0.216 
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Clearly the best overall fit is the simultaneous fit with two a poles a and a', Solution 
A. Its x 2 /dpt is lower than x 2 /dpt for separate fits, and much lower than the y^/dpt for 
simultaneous fits with a single a pole. From this we can conclude that the hypothesis that 
the width of a depends on nucleon spin has been invalidated. Equally invalidated is the 
hypothesis of a single a pole contributing to both \S\ 2 H and \S\ 2 H amplitudes. Since the 
Solution B has a larger x 2 /dpt close to x 2 jdpt for separate fits, we propose to reject it as 
unphysical. As the result of this x 2 criterion we are left with two a states a and a 1 from the 
Solution A. A possible physical interpretation of these states will be discussed in Section VI. 

C. Spin dependence of free parameters. 

It is of interest to examine the dependence of free parameters on nucleon spin. In this 
discussion we will restrict attention mostly to the Solution A. 

A remarkable feature of all fits is that in all amplitudes the contribution of the higher 
mass and broader a dominates the contribution of a 1 . In Table IV we present the real and 
imaginary parts Di and D 2 of the coupling of a' Breit-Wigner amplitude. We see that for 
recoil transversity "up" amplitudes 1 and 2, the real part D\ is negative and imaginary part 
D 2 is positive. The opposite is true for D 1 and D 2 for the recoil nucleon transversity "down" 
amplitudes 1 and 2. The average absolute magnitude \D\ av also depends on recoil nucleon 
transversity: 



Solution A 0.08 0.40 0.62 0.60 

Solution B 0.18 0.35 0.61 0.61 

We see that the coupling of a' is weakest in the amplitude IS^E, Solution 1. It is strongest 
in spin down amplitude \S\ 2 H in both solutions. 

In Table V we show the real and imaginary parts C\ and C 2 of the coupling of / (980) 
together with the fitted values of mj and Tf. We see that C\ and C 2 are both positive in 
both solutions for spin up amplitude l^ 2 !!. In contrast, C\ is negative in both solutions 
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for spin down amplitude [S'pS. The fitted values of m/ and Tf are in agreement with 
expectations for / (980). 

Finally, in Table VII we present the real and imaginary parts B\ and B 2 of coherent 
background and the overall normalization constant. B\ and B2 are both positive for both 
solutions of spin up amplitude l^pE. For spin down amplitude \S\ 2 Ti the real part B\ is 
negative while the imaginary part B 2 remains positive. The values of normalization constant 
Ns show that a couples most strongly to Solution 1 of the spin up amplitude |5| 2 £. 

VI. POSSIBLE CONSTITUENT INTERPRETATIN OF THE SCALAR STATES 

<r'(670) AND <r(786). 

In the usual quark model meson resonances are qq states. The mass M of the qq states 
increases with its angular momentum L as M = M (2n + L) where n is the degree of radial 
excitation. The lowest mass scalar mesons are 3 Pq states with masses expected to be around 
1000 MeV or higher. The masses of <r'(670) and <r(786) are too low to be qq states. 

The fact that cr'(670) and <r(786) are not observed in reactions 77 — > and 77 — > 

7r°7r° suggests a gluonium interpretation of these states. Since gluons do not couple directly 
to photons, we expect the a' and a states to be absent in 77 — > tttt reactions. This conclusion 
is supported by PLUTO and DELCO data |pO| , |3~Tf and more recent results on 77 — ► n - -" 



7T 



Ellis and Lanik discussed the couplings of scalar gluonium a on the basis of the low-energy 
theorems of broken chiral symmetry and scale invariance, implemented using a phenomeno- 
logical Lagrangian |l9[. They obtained for a — ► 7r + 7r _ decay the following partial width 



r( a - 7T+7T-) = (6.1) 
4o7rCrn 



where Go =< ^{(o^s/^F^F^lO > is the gluon condensate term J33] parametrizing the 
nonperturbative effects of QCD. The numerical values were estimated by the ITEP group 
to be G ~ 0.012 (GeV) 4 or up to G « 0.030 (GeV) 4 in later calculations mjS5 



Several recent estimates of G all agree on the values around G ~ 0.020 (GeV) 4 |]36|-p8 
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It is very interesting to note that when we take Go = 0.015 (GeV) 4 , the Ellis-Lanik 
theorem (6.1) predicts T(a — > 7r + 7r _ ) = 133 MeV for the a mass m a = 786 MeV and 
r(cr' — > 7T + 7r _ ) = 60 MeV for the a' mass m a i = 670 MeV. This result is in perfect agreement 
with Eq. (5.12) where T a = 130 ±47 MeV and IV = 59 ±58 MeV. However, T a and T a , are 
the full widths of <r(786) and <r'(670), so these results on T(cr — > 7r + 7r~) and T(a' — > 7r + 7r~) 
represent upper limits. When we use Go — 0.020 (GeV) 4 , which is the average of the 



latest values for Go [pq-pql, we get for the partial width T(a — > 7r + 7r~) = 0.75r CT and 
r(cr / — > 7r + 7r _ ) = 0.75FO-/ with a very reasonable branching fraction of 75% for the 7r + 7r~ 
channel. From this agreement with Ellis-Lanik theorem we can conclude that a' (670) and 
cr(786) are both best understood as the lowest mass scalar gluonium states ++ (gg)- 

For further interpretation of the cr'(670) and <r(786) states we turn to Bjorken's classifi- 
cation of gluonium states with no spatial excitation p2"|-f24}j. In our discussion we will follow 



Ref. [25 



Consider gluon field strength tensor 

= - d u Al + g s f aPl A^Al (6.2) 

where A° are QCD gluon gauge fields. The tensor F" transforms like an octet under SU(3) C 
colour symmetry group. In analogy to electrodynamics we define colour electric fields Ef 
and colour magnetic fields Bf. 



The colour electric fields correspond to J PC = 1 states while the colour magnetic fields 
correspond to J PC = l + ~ states p5| . Since gluons are Bose particles, their total wave 
function must be symmetric with respect to the interchange of all labels. We restrict our 
consideration to s-waves with no spatial excitations which have a symmetric space wave 
function. Thus the product of the spin and colour wave functions must be also symmetric. 
Since bound gluons are off-shell, we must also admit the longitudinal spin component. Then 
there are two distinct modes of lowest mass scalar gluonium ++ (gg) P5| : colour electric 
mode gluonium 
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8 3 

\0 ++ EE>= £ £ ScpSijlEfEf > (6.4) 



a, /3=1 i,j=l 

and colour magnetic mode gluonium 

8 3 



\0 ++ MM>= £ EWilW> (6-5) 

a,/3=l i,j=l 

In (6.4) and (6.5) the indices a, (3 = 1, . . . ,8 stand for colour spin, while the indices i,j = 
1, . . . , 3 stand for the polarization of the field which is related to the third component of the 
spin of the corresponding off-shell gluon |25| . 

It is expected that the magnetic mode (6.5) is heavier than the electric mode (6.4) 
22] , |23|p5| . Magnetic mode scalar gluonium is also heavier than electric mode scalar gluonium 



in MIT bag model ||39|| . On the basis of these considerations we propose to identify the lower 
mass er'(670) state with the colour electric mode \0 ++ EE >, and the higher mass cr(786) state 
with the colour magnetic mode \0 ++ MM > of the lowest mass scalar gluonium ++ (gg). 



VII. SUMMARY. 



We have performed a new model independent amplitude analysis of reaction 7i~p — > 
7T~7T + n at 17.2 GeV/c extending the range of dipion mass from 600-880 MeV of the earlier 
analysis Jl7 , 18| to 580-1080 MeV in the present analysis. In this study we focused on 



the resonant structure of the S- wave amplitudes |5| 2 S and \S\ 2 H in this new mass range. 
Our purpose was to understand the role of / (980) interference and to determine better the 
resonance parameters of <r(750) resonance from fits to both transversity amplitudes \S\ 2 H and 
l^pE. This work is the first attempt to determine resonance parameters of a resonance from 
fits to production amplitudes of opposite nucleon transversity. In the absence of previous 
experience, we used several fitting procedures to understand the resonant structure of the 
S-wave amplitudes. 

The amplitude \S\ 2 Y> (recoil nucleon transversity "up") shows a clear resonant peak at 
750 MeV below 880 MeV, an enhancement above 900 MeV followed by a rapid descent to a 
dip at 1010 MeV. The enhancement is more pronounced in the Solution 2. The enhancement 
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and the dip are explained by the interference of /o(980) with coherent background and the 
contribution from a(750) resonance. The interference of /o(980) and cr(750) does not lead 
to a high mass and a broad width of a as suggested in Ref . pO| , pT]| . 

The amplitude |S| 2 E (recoil nucleon transversity "down") shows a broad structure and 
a dip at 1010 MeV. The broad structure is a composite of a contribution from cr(750) 
resonance and enhancement due to interference of /o(980) with background and o. The 
broad structure arises as the two contributions have comparable magnitudes. The resonance 
behaviour of \S\ £ was not recognized in the earlier analyses [|T7l , |T8l because the mass range 
was limited only to 600-880 MeV. As the result the amplitude |S| 2 £ was described in |T7|JT8[ 
as "nonresonating" . 

The above description of the structure of S-w&ve amplitudes was first reached by using 
separate fits to the amplitudes \S\ 2 Y, and j^^E. However the separate fits suffer from the 
difficulty that the width of a is narrow 110 or 76 MeV in Solutions 1 and 2 of amplitude 
|5| 2 E but it is broad at 217 and 233 MeV in Solutions 1 and 2 of amplitude jS'PE. 

Next we tried to ascertain whether this apparent dependence of a width T a on nucleon 
spin is a real physical effect or an artifact of separate fits to |5| 2 E and (S'pE. To this end 
we have first performed a simultaneous fit to \S\ 2 H and \S\ 2 Ti with a common a pole. The 
differences in the shape of mass distributions |S"| 2 E and \S\ 2 H were expected to arise purely 
from the coherent background. The result of the best fit gives 

m CT = 775 ± 17 MeV , T a = 147 ± 33 MeV (7.1) 



We note that the previous determinization |T8| of m a and IV for |,S| £ alone and without 



/o(980) interference gave somewhat smaller mass and a narrower width. 

m a = 753 ± 19 MeV , r CT = 108 ± 53 MeV (7.2) 

The values of x 2 / dpt in the best simultaneous fit show a good fit but nevertheless they are 
higher than the x 2 '/dpt from the best separate fits. Thus the possibility of the dependence 
of Y a on nucleon spin has not been eliminated by these fits. Consequently we turned to 
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examine the hypothesis that the width IV does depend on nucleon spin in a new set of 
simultaneous fits to \S\ 2 H and IS'pE. If we denote a the narrow a contributing to \S\ 2 T, and 
a the broader a contributing to 1 5"| 2 S, the self-consistency requires that some a contributes 
to \S\ 2 H and some a contributes to l-S'pE. So instead of one a pole with a width dependent 
on nucleon spin, we were now looking for two a poles of similar masses with different widths 
and different relative contributions to the amplitudes \S\ 2 H and \S\ 2 T* to explain the results 
of separate fits. 

The results of the simultaneous fits with two a poles turned out differently. We did 
find two a poles with different widths, but they also have considerably different masses and 
did not distinguish themselves by contributing preferentially to one or the other S'-wave 
amplitude. We called a the a pole with higher mass and broader width, and a' the a pole 
with lower mass and narrower width. We found that a is the dominating contribution in 
both solutions in both amplitudes \S\ 2 H and \S\ 2 T,. The contribution of a' is weaker in all 
amplitudes, but it does depend on the nucleon spin. It is weak in the spin "up" amplitudes 
l-S'pE, and much stronger in the spin "down" amplitudes j^pE. 

We found two solutions for the two poles a and a'. The Solution A is much prefered on 
the basis of very low x 2 /dpt. The averaged resonance parameters are 

m a = 786 ± 24 MeV , T a = 130 ± 47 MeV (7.3) 

77V = 670 ± 30 MeV , T a , = 59 ± 58 MeV 

The Solution A for a and a' has x 2 /dpt lower than x 2 /dpt for separate fits or simultaneous 
fits with single a pole. On the basis of this x 2 criterion we recommend that the Solution A 
for a and o' poles be accepted as the actual physical explanation of the resonant structure 
of amplitudes |5| 2 E and \S\ 2 Y,. 

In Section VI we showed that both a and a' poles meet several criteria for gluonium 
candidates: their masses are too low to be a qq scalar states, they are not observed in 
77 — > mi reactions, and they both satisfy Ellis-Lanik theorem (6.1) with the same value of 
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gluon codensate Go- On the basis of Bjorken's classification of gluonium states p2|-f25 



we 



proposed to identify cr'(670) and cr(786) with the colour electric and colour magnetic modes 
of lowest mass scalar gluonium ++ (gg). 

The CERN measurements of pion production on polarized targets have opened a whole 
new approach to experimental hadron spectroscopy by making possible the study of res- 
onance production on the level of spin amplitudes. Our results emphasize the need for a 
dedicated and systematic study of various production processes on the level of spin ampli- 
tudes measured in experiments with polarized targets. 
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FIGURES 

FIG. 1. Mass dependence of physical solutions for moduli squared of S- wave and P-wave nu- 
cleon transversity amplitudes and cosines of their relative phases in reaction ir~p — ► 7r~7r + n at 
17.2 GeV/c and momentum transfers —t = 0.995 — 0.2 (GeV/c) 2 . The results are in the t-channel 
dipion helicity frame. 

FIG. 2. Mass dependence of unnormalized amplitudes \S\ 2 T, and \S\ 2 H measured in 
7r-p T -► Tr-7r + n at 17.2 GeV/c and -t = 0.005 - 0.20 (GeV/c) 2 . Figures (a) and (b) show 
the results using the Monte Carlo method for amplitude analysis of Data Set 1 (Ref. [17]) and 
Data Set 3 (this paper). Figures (c) and (d) show the results using the x 2 minimization method 
for amplitude analysis of Data Set 1 (Ref. [5,18]) and Data Set 3 (Ref. [10]), respectively. 

FIG. 3. Four solutions for the S"-wave partial-wave intensity Is in the reaction ir~p — ► 7r~7r + n 
at 17.2 GeV/c and -t = 0.005 - 0.20 (GeV/c) 2 . 

FIG. 4. The fits to amplitudes |5| 2 S and |£| 2 E. (a) The separate fits with separate single 
a pole using the parametrization (3.9). (b) The simultaneous fit with a common single a pole 
using the parametrization (3.9). Figures (c) and (d) show Solution A and B, respectively, of the 
simultaneous fit with two common a poles using the parametrization (5.9) and (5.10). The Solution 
A in (c) gives the best \ 2 /dpt. 
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TABLES 

TABLE I. Mass and width of cr(750) and /o(980) from separate fits to solutions 1 and 2 of 
amplitude |S*| 2 S, and to solutions 1 and 2 of amplitude |5| 2 E. 



Solution 


m a 


L(j 


m f 




X 2 /dpt 




(MeV) 


(MeV) 


(MeV) 


(MeV) 




T 


773 ± 12 


110 ± 25 


956 ± 17 


81 ± 40 


0.206 


2 


775 ± 16 


76 ± 37 


973 ± 17 


57 ± 25 


0.240 


1 


788 ± 33 


217 ± 53 


982 ± 10 


42 ± 39 


0.294 


2 


807 ± 45 


233 ± 109 


988 ± 10 


16 ± 25 


0.148 



TABLE II. Mass and width of cr(750) and / (980) from simultaneous fits to |Sf£ and |5| 2 S 
with a common single a pole. The combination of solutions of |S*| 2 S and |S| 2 £ is given in the 
parentheses on the left with the notation as in Table I. The x 2 /dpt shown is the average of x 2 /dpt 
for |5| 2 £ and |,S| 2 £ solutions. 



Fit 


m a 


Lcr 


m f 




X 2 /dpt 




(MeV) 


(MeV) 


(MeV) 


(MeV) 




(11) 


771 ± 13 


161 ± 22 


980 ± 10 


49 ± 21 


0.386 


(12) 


771 ± 14 


145 ± 25 


970 ± 19 


60 ± 28 


0.292 


(21) 


780 ± 21 


156 ± 43 


980 ± 9 


47 ± 25 


0.370 


(22) 


780 ± 18 


125 ± 40 


975 ± 45 


20 ± 20 


0.272 
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TABLE III. Mass and width of two a poles from simultaneous fits to \S\ 2 H and \S\ 2 H using 
parametrization (5.9), (5.10) for |S*| 2 S and corresponding parametrization for |5| 2 S. Values with- 
out errors were obtained in a constrained fit. The x 2 /dpt shown is the average of \ 2 I dpt for |S*| 2 S 
and \S\ 2 H solutions. 



Amplitudes 


m a 


r CT 


m c i 


rv 


X /dpt 




(MeV) 


(MeV) 


(MeV) 


(MeV) 




Solution A 












(11) 


775 


118 


676 


73 


0.274 


(12) 


775 ± 19 


118 ± 34 


676 ± 44 


73 ± 71 


0.168 


(21) 


I I'd ± 2 ( 


\X>( ± DO 


066 ± 12 


36 ± 38 


n 1 to 
U.LlZ 


(22) 


795 ± 29 


120 ± 57 


665 ± 18 


55 ± 51 


0.110 


Solution B 












(11) 


761 ± 27 


140 ± 69 


717 ± 43 


72 ± 79 


0.268 


(12) 


770 ± 23 


123 ± 48 


709 ± 44 


69 ± 77 


0.168 


(21) 


770 


123 


709 


69 


0.246 


(22) 


770 


123 


709 


69 


0.180 
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TABLE IV. Spin dependence of the a' coupling constants D\ and D 2 denned in (5.5). D\ and 
D2 are understood for amplitudes 1 and 2. 



Amplitude 


m a ,T a 


(MeV) 


m a : , T a i 


(MeV) 




D 2 


T(T1) 


(775, 


118) 


(676, 


73) 


-0.06 ±0.10 


0.05 ±0.10 


1(12) 


(775, 


118) 


(676, 


73) 


-0.06 ±0.20 


0.06 ±0.15 


2(21) 


(799, 


167) 


(666, 


36) 


-0.36 ± 0.36 


0.20 ±0.25 


2(22) 


(795, 


120) 


(665, 


55) 


-0.31 ±0.36 


0.25 ± 0.29 



1(11) 

1(21) 
2(12) 
2(22) 



(775, 118) 
(799, 167) 
(775, 118) 
(795, 120) 



(676, 73) 
(666, 36) 
(676, 73) 
(665, 55) 



0.80 ± 0.60 
-0.06 ±0.33 
0.65 ±0.88 
0.10 ±0.52 



-0.38 ±0.53 
-0.34 ±0.36 
-0.27 ±1.37 
-0.47 ±0.56 



TABLE V. Spin dependence of the /o(980) coupling constants C\ and C 2 defined in (5.5) and 
the fitted values for m/ and Tf. C\ and C 2 are understood for amplitudes 1 and 2. The ordering 
of masses is as in Table IV. 



Amplitude d C 2 m f (MeV) T f (MeV) 

1(11) 0.15 ±0.28 0.65 ±0.20 968 ± 16 95 ± 29 

T(T2) 0.30 ±0.55 0.60 ±0.26 958 ± 18 91 ±40 

2(21) 0.92 ± 1.04 1.34 ± 1.44 981 ± 7 35 ± 24 

2(22) 1.38 ± 1.05 0.89 ± 1.41 973 ± 16 58 ± 26 

l(Tl) -1.37 ±0.71 -0.21 ±0.86 968 ± 16 95 ± 29 

1(21) -0.83 ±0.45 0.35 ±0.56 981 ± 7 34 ± 24 

2(12) -1.09 ± 1.29 0.33 ± 1.30 958 ± 19 91 ±41 

2(22) -0.66 ±0.85 0.56 ± 0.71 974 ± 15 58 ± 26 
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TABLE VI. Spin dependence of coherent background B\ and B 2 and the normalization con- 
stant Ns defined in (5.5), (5.9) and (5.10). B±, B2 and N$ are understood for amplitudes 1 and 
2. The ordering of masses as in Table IV. 



Amplitude B x B 2 N s 

1(11) 0.82 ±0.19 0.27 ±0.13 2.81 ± 0.74 

T(T2) 0.81 ±0.29 0.22 ±0.35 2.91 ± 0.91 

2(21) 1.59 ±0.64 0.36 ±0.75 1.42 ±0.82 

2(22) 1.82 ±0.73 0.15 ± 1.26 1.26 ±0.82 

1(11) -0.92 ±0.60 0.44 ±0.42 0.76 ± 0.61 

1(21) -0.12 ±0.25 0.35 ±0.44 1.92 ± 1.56 

2(T2) -1.10 ± 1.42 0.61 ±0.98 0.80 ± 1.26 

2(22) -0.47 ±0.74 0.89 ± 0.65 1.15 ± 1.06 
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